TUM-HEP-511/03 
HD-THEP-03-26 



e'/e at the NLO: 10 Years Later 



Andrzej J. Buras a and Matthias Jamin 6 * 

, a Physik Department, Technische Universitiit Miinchen 

(N . D-85748 Garching, Germany 

b Institut fur Theoretische Physik, Universitat Heidelberg, 

£NJ ' Philosophenweg 16, D-69120 Heidelberg, Germany 

(N ' 



(N 
> 



Abstract 

During the last four years several parameters relevant for the analysis of the CP-violating 
ratio e'/e improved and/or changed significantly. In particular, the experimental value of e'/e 

o 

and the strange quark mass decreased, the uncertainty in the CKM factor has been reduced, 
and for a value of the hadronic matrix element of the dominant electroweak penguin operator 
Q%, some consensus has been reached among several theory groups. In view of this situation, 
Q_i' ten years after the first analyses of e'/e at the next-to-leading order, we reconsider the analysis 

of e je within the SM and investigate what can be said about the hadronic Qq matrix element 
of the dominant QCD penguin operator on the basis of the present experimental value for e'/e 
^ ■ and todays values of all other parameters. 

Employing a conservative range for the reduced electroweak penguin matrix element R% = 
1.0 ±0.2 from lattice QCD, and present values for all other input parameters, on the basis of the 
current world average for e'/e, we obtain the reduced hadronic matrix element of the dominant 
QCD penguin operator R 6 = 1.23 ± 0.16 implying (<2 6 )o DR ( m c) « - 0.8 (Q 8 )^ DR (m c ). We 
compare these results with those obtained in large-iV c approaches in which generally Rq ~ Rs 
and (Q6)o DR ( m c) is chirally suppressed relatively to (Qg)^ DR (m c ). We present the correlation 
between Rq and Rg that is implied by the data on e'/e provided new physics contributions to 
e'/e can be neglected. 
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1 Introduction 



The ratio e'/e that parametrises the size of direct CP violation with respect to the indirect CP 
violation in Kl — ► vr7r decays has been the subject of very intensive experimental and theoretical 
studies in the last three decades. After tremendous efforts, on the experimental side the world 
average based on the recent results from NA48 [1,2] and KTeV [3,4], and previous results from 
NA31 [5,6] and E731 [7], reads 

e'/e = (16.6 ± 1.6) • 1(T 4 (2003) . (1.1) 

A recent discussion of the experimental issues involved in these analyses can be found in ref. [8]. 

On the other hand, the theoretical estimates of this ratio are subject to very large hadronic 
uncertainties. While several analyses of recent years within the Standard Model (SM) find results 
that are compatible with (1.1) [9-18], it is fair to say that even ten years after the first analyses 
of e'/e at the next-to-leading order (NLO) [19,20], the chapter on the theoretical calculations 
of e'/e is certainly still open. A full historical account of the theoretical efforts before 1998 can 
for example be found in the reviews [21,22]. 

It should be emphasised that all existing analyses of e'/e use the NLO Wilson coefficients 
calculated by the Munich and Rome groups in 1993 [19,20,23-26], but the hadronic matrix 
elements, the main theoretical uncertainty in e'/e, vary from paper to paper. Nevertheless, apart 
from the hadronic matrix element of the dominant QCD penguin operator Qq, in the last years 
progress has been made with the determination of all other relevant parameters, which enter the 
theoretical prediction of e'/e. In view of this situation, we think it is legitimate to reconsider 
the analysis of e'/e within the SM and to investigate what can be said about the hadronic Qq 
matrix element on the basis of the present experimental value of e'/e and todays values of all 
other parameters. In doing this we assume that any potential new physics contributions to e'/e 
can be neglected. We will return to this point at the end of our paper. 

The phenomenological approach towards the values of the hadronic matrix elements pre- 
sented in this work is a continuation of our approach in [19], where we extracted the matrix 
elements of current-current operators Qi (i = 1 — 4) from the CP conserving data. In that case 
the assumption about the absence of significant new physics contributions to the relevant am- 
plitudes was not necessary. On the contrary, the present analysis has this particular limitation. 
Nevertheless, as we shall see below some insight into the pattern of the size of the hadronic 
matrix elements of penguin operators can be gained. 

In order to set the scene for our investigation let us recall that after NA31 in 1993 [6] and 
KTeV in 1999 [4] presented surprisingly high values for e'/e, the 1999 world average including 
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the E731 value [7] was 

e'/e = (21.8 ± 3.0) • 10~ 4 



(1999) , 



(1.2) 



being roughly 30% higher than the 2003 value. Within our 1999 analysis of e'/e [17] that was 
performed using matrix elements of penguin operators Qq and Qg in the ballpark of the values 
obtained in the large- N c approach of [27-29], it was essentially impossible to reproduce the 
value in (1.2). We found e'/e typically a factor of two to three below the 1999 experimental 
value. 

As we shall see below, the confrontation of the SM with the 2003 experimental value for e'/e 
as given in (1.1), does not resemble the 1999 situation presented by us in [17]. The following 
important facts are responsible for this change: 

• The experimental value of e'/e decreased approximately by a factor of 1.3. 

• The value of m 8 , relevant for e'/e within the approach of [27-29], decreased by roughly 
15%, enhancing the theoretical value of e'/e by roughly a factor of 1.3. 

• The isospin breaking parameter Otb decreased from roughly 0.25 to 0.06, enhancing e'/e 
by another factor of 1.4. 

• The value for the hadronic matrix element of the dominant electroweak penguin operator 
(Qs) became more robust due to general consistency of results from lattice QCD and other 
analytic non-perturbative approaches to be discussed below. 

• Finally, the value of the CKM factor lm.X t became more precise, although it did not change 
significantly compared to our 1999 value. 

As a result of all these changes, our SM value for e'/e presented below turns out to be fully 
consistent with the experimental value given in (1.1) for values of the hadronic matrix elements 
(Qq) and (Qs) which do not differ by more than 20% from the ones obtained within the large-A^ c 
approach of [27-29]. On the other hand, the chiral suppression of (Qe)o with respect to (Qs)2 
that is characteristic for the latter approach, is only partially supported by the data. 

In the next section, we recollect the basic equations which are required for the analysis of 
e'/e within the SM. In section 3, we present the details and numerical results of this work, in 
section 4 we compare our results with those obtained in various approaches, and in section 5 we 
summarise and draw some conclusions. 
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2 Basic Formulae 



The central formula for e'/e of refs. [17,19,30] takes the following form: 



— = ImAt • F £ >(x t ) 

e 



where 



x t = V t *V td , x t 



m t 



and 

F £ ,{x t ) = P + Px X Q {x t ) + Py Yo(x t ) + P z Z (x t ) + P E E (x t ) . 
The functions Xq, Yq, Zq and Eq are given by (m t = 165 GeV) 

= 1.505 , 



X (xt) 
Y (xt) -- 



x t 



x t + 2 3x t - 6 

+ ~ — l n Xt 



x t -l (x t - 1) 



xt-4 



+ 



3x* 



In x t 



0.962 , 



1 



Z (x t ) = -- kix t + 



[xt-l ' {xt-lf 

18x 4 - 163xf + 259x 2 - 108x t 



+ 



144(x t - l) 3 
32x 4 - 38x 3 - 15x 2 + 18x t 



+ 



72(x t - l) 4 
x 2 

E (x t ) = -^lnx t + ' 



lnxt = 0.664 



2 , x 2 (15 - 16x t + 4x t 2 ) x t (18 - llx t - x\) 

In Xt -\- 



6(1 -x t ) 4 



12(1 -x t ) 3 



= 0.271 . 
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-3.574 
0.574 
0.403 
0.714 
0.213 


16.552 
0.030 
0.119 
-0.023 
-1.909 


1.805 



-12.510 
0.550 


-3.602 
0.564 
0.392 
0.766 
0.202 


17.887 
0.033 
0.127 
-0.024 
-2.017 


1.677 



-13.158 
0.589 


-3.629 
0.554 
0.382 
0.822 
0.190 


19.346 
0.036 
0.134 
-0.026 
-2.131 


1.538 



-13.855 
0.631 



Table 1: The coefficients rf \ rf^ and rf of formula (2.8) for various A^L in the NDR scheme. 



,(8) 



MS 
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The coefficients Pi are given in terms of the non-perturbative parameters Rq and Rs as 
follows: 



rf +rfi? 6 + rfi?g, (2.8) 



where the coefficients rf 1 ^ , rf^ and rf^ comprise information on the Wilson-coefficient functions 
of the AS = 1 weak effective Hamiltonian at the next-to-leading order [23-26], and their nu- 
merical values for different values of at a = m c in the NDR renormalisation scheme are 

MS ^ 

displayed in table 1. This table updates the corresponding ones presented in [17,30]. 

(1/2) (3/2) 

In terms of the hadronic S-parameters Bq and B% ' , as well as the strange quark mass, 
the non-perturbative parameters Rq and Rs are defined as 



r> - R (l/2) 



121 MeV 



_m s {m c ) + m d {m c ) 



Rs = b?i*> 



121 MeV 



m s {m c ) + m d {m c ) 



(2.9) 



In the large- A^ approach of [27-29], the matrix elements of the dominant QCD-penguin (Qq) 
and the dominant electroweak penguin (Qs) operator, are then given by [19,31] 



Wo 



vf {Fk ~ F ^ (i^rlev) 2 ^ = - °- 597 • R « GeV3 ' (2 - 10) 



(Q 8 ) 2 = ^f J™^\ R 8 = 0.948 • R 8 GeV 3 , (2.11) 

where for {Qs)2, w e have neglected a small 1% correction which is independent of m s , and the 
numerical values have been obtained by employing values for F n , Fk and ijik as given in the 
current PDG [32]. Let us remark that the overall normalisation of our matrix elements in (2.10) 
and (2.11) differs by a factor \/3/2 from another convention used in the literature. This factor, 
however, compensates with a corresponding factor in our definition of the isospin amplitudes to 
provide the same physical results. 

(1/2) (3/2) 

In the strict large-A^ c limit, B$ ' = Bg — 1, the /i-dependence of the matrix elements 
{Qe)o an d (Qs}2 is fully governed by the /i-dependence of the quark masses m S; d(/x), and, as 
the mixing with other operators can be neglected in this limit, the evolution of the Wilson 
coefficients Cq(/j.) and Cg(/x) precisely cancels the ^-dependence of (Qe)o and (Q$)2 such that 
e'/s is /U-independent as it should be. As our numerical analysis of ref. [19] demonstrated, B^ 2 ^ 

(3/2) 

and B 8 remain independent of \x to a very good accuracy also outside the strict large- N c limit, 
provided these parameters do not deviate considerably from unity. 

The formulae (2.9), together with the strict large-A" c limit, also imply Rq = Rg, so that in 
this limit there is a one-to-one correspondence between e'/e and Rq = Rs for fixed values of the 
remaining parameters. Consequently, in this case, only for certain values of m s (m c ) is one able 
to obtain the experimental value for e'/e. This has been emphasised in particular in [33], and 
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we shall return to this point in the next section. Thus in this approach a concrete prediction 
for e'/e can be made provided the value m s (m c ) has been determined by some other method. 
Moreover, once m s (m c ) is known, also Rq, R$, {Qe)o and (Qs}2 are known, but they always 
satisfy the relations: 

f = l, ^ = -0.63. (2.12) 

Rs \V8/2 

We shall compare the relations (2.12) with our results for (Qq)o and (Qs)2 obtained by various 
methods below. 



3 Numerical Analysis 
3.1 Preliminaries 

Before we come to a detailed presentation of our numerical analysis, let us first briefly discuss 
the origin of the values for the dominant input parameters that we have used in our analysis. 

The dependence of e'/e on the quark mixing or CKM matrix elements appears through the 
parameter ImA^. Fits determining the elements of the CKM matrix have been the subject of 
extensive efforts in the last years. The present status has been summarised very recently in 
the review [34]. Employing the best fit incorporating all known constraints within the SM, the 
resulting value for ImAt is then found to be [35]: 

ImA t = (1.31 ±0.10) • 1(T 4 . (3.1) 

The next parameter which deserves an explicit mentioning since e'/e is rather sensitive to 
its precise value, is the quantity SIib, parametrising isospin-breaking corrections which result 
from strong isospin violation (m u ^ m^) and electromagnetic corrections. For this parameter, 
we employ the result of the very recent analysis [36] (see also refs. [37,38]): f^iB = 0.06 ± 0.08. 
The shift from S7ib = 0.25, that was used in our previous analysis of e'/e [17], to the new result, 
has the largest impact for the theoretical value of e'/e. 

An important fundamental QCD input parameter is the value of the strong coupling constant 
a s . As pointed out in refs. [17,39], e'/e is roughly proportional to the QCD parameter Aj^g in a 
four-flavour theory, when this parameter is 0(350 MeV). For the strong coupling, we have used 
the recent average a s {M z ) = 0.119 ± 0.002, corresponding to A^ = 340 ± 30 MeV. 

Another parameter is the top quark mass that enters the Inami-Lim functions given above 
and the fit of ImAt. However, wit is already well known and its precise value is less important. 
On the other hand, the value of the strange quark mass m s plays an important role in the 
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large-iV c approach discussed in the previous section. For the quark masses, we shall then use 
the following values: 

m t {m t ) = 165 ± 5 GeV , m s {m c ) = 115 ± 20 MeV . (3.2) 

The top mass evaluated at the scale m% results from converting the corresponding pole mass 
value [32] , and the strange quark mass is an average over recent determinations of this quantity 
[40-44]. Due to our method of extracting the matrix elements of current-current operators 
Qi (i = 1 — 4) from the CP conserving data [19], the strange mass is required at the scale 
m c = 1.3 GeV and our central value would correspond to m s (2 GeV) = 100 MeV. 

As stressed by several authors in the literature, the strong m s dependence of {Qe)o and (Qs)2 
given in the large--/V c approach is somewhat artificial. In fact, the lattice studies [45] and the 
discussion in refs. [9,13,46] indicate that these matrix elements are only weakly dependent on the 
actual value of m s since they are related to the quark condensate via the GMOR relation [47, 48] . 
Here we shall follow first a different approach. We will directly deal with the parameters Rq and 
R$ of eq. (2.9), and only at the end of our numerical analysis, we shall discuss the implications 
for and i?^. 

3.2 Present knowledge of R 8 

As can be inferred from the formulae presented in the last section, in order to be able to calculate 
s' /e, we still require some input for Rq and Rg which parametrise the dominant contributions 
from QCD and electroweak penguin operators respectively. Concerning the parameter R$, in 
the last years progress has been achieved both in the framework of lattice QCD [49-51] as well 
as with analytic methods to tackle QCD non-perturbatively [9,52-55]. 

The current situation of the status of {Qs)2 calculations has been summarised nicely in the 
recent article [55]. At present, the most precise determination of (Qs)2 comes from the lattice 
QCD measurement [51], corresponding to Rg = l.OrtO.l. The lattice studies of refs. [49,50] find 
results in complete agreement, however, with somewhat larger uncertainties. The analytic non- 
perturbative approaches of refs. [52-55], employing either the large- N c expansion or dispersion 
relations, on the other hand, find results that are generally larger than the lattice findings, 
although certainly compatible within the uncertainties. However, these analyses have been 
performed in the chiral limit, and as was demonstrated in [55], chiral corrections tend to decrease 
Rg. Therefore, in order to be more conservative, and in view of the fact that there could be 
additional effects of chiral logarithms when relating the matrix elements computed on the lattice 
to the physically relevant ones, in what follows we shall use: 

R$ = 1.0 ± 0.2 , {Qs)2 DR (. m c) = (0.95 ± 0.19) GeV 3 . (3.3) 



7 



Rg 


Rq 


e'/e 


ImAj 






m t 


m s 


total 


0.8 


1.143 


+0.070 
-0.070 


-0.051 
+0.060 


-0.061 
+0.065 


+0.092 
-0.078 


+0.018 
-0.017 


+0.007 
-0.012 


+0.147 
-0.133 


1.0 


1.222 


+0.070 
-0.070 


-0.051 
+0.060 


-0.061 
+0.064 


+0.099 
-0.084 


+0.024 
-0.023 


+0.007 
-0.012 


+0.152 
-0.138 


1.2 


1.301 


+0.070 
-0.070 


-0.051 
+0.060 


-0.061 
+0.064 


+0.107 
-0.090 


+0.029 
-0.028 


+0.007 
-0.012 


+0.158 
-0.142 


2.2 


1.695 


+0.070 
-0.070 


-0.051 
+0.060 


-0.059 
+0.062 


+0.143 
-0.121 


+0.058 
-0.057 


+0.007 
-0.012 


+0.190 
-0.170 



Table 2: Rq as a function of Rs for four relevant values. In addition, we give the variation of 
Rq while varying all other input parameters. For a more detailed explanation see the text. 

3.3 Determining R 6 

The situation is less clear concerning the QCD penguin contribution parametrised by Rq. There- 
fore, we postpone a discussion of the present status of Rq, and we shall take a different approach. 
Putting together all the information presented thus far, let us investigate what can be said about 
our expectations for Rq. A summary of this analysis is presented in table 2. The second column 
shows central values for Rq corresponding to the four inputs for Rg as given in the first column. 
The remaining columns display the variation of Rq while varying the dominant input parame- 
ters, and the last column contains the total uncertainty when all individual errors are added in 
quadrature. 

A graphical representation of our results of table 2 in the spirit of the "e'/e-path" of ref. [56] 
(see also refs. [9,46]) is given in figure 1. The solid straight line corresponds to the central values 
given in table 2, whereas the short-dashed lines are the uncertainties due to a variation of the 
input parameters. (Let us note that all our errors should be considered as la deviations.) Next, 
the vertical long-dashed lines indicate the lattice range for Rg [51], whereas the ellipse describes 
the correlation between Rq and Rg implied by the data on e'/e when taking into account the 
more conservative constraint on Rg given in eq. (3.3). From the allowed region within the ellipse, 
we can infer our central result for Rq, which reads: 

Rq = 1.23 ± 0.16 , (<2 6 >o DR (™c) = - (0.73 ± 0.10) GeV 3 . (3.4) 

Finally, the full triangle, circle and diamond in figure 1 represent the central results of refs. [13, 
54], [10] as well as [9] respectively, which will be discussed in more detail below, and the dashed- 
dotted line shows the strict large- iV c relation Rq = Rg. 

Comparing (3.4) with (3.3), we observe that if the matrix element of Qg is indeed given by 
(3.3) than the data on e'/e seem to require 

^ = 1.23 ±0.29, <Q 6 >r R (™ C ) = - (0.77 ± 0.19) (Qg)^ DR (m c ) . (3.5) 
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Figure 1: fig as a function of R%. For a detailed explanation see the text. 

Due to substantial uncertainties, these results are compatible with the large-iV c expectations 
in (2.12) but the strong chiral suppression of the matrix element (Q6)q DR (iti c ) with respect to 
(Q8)2 DR ( m c) found in the latter approach, although visible in (3.5), is not fully supported by 
the data. However, we would like to emphasise that this result is a consequence of (3.3). We 
will return to this point in the next section, where we shall compare the results in (3.3), (3.4) 
and (3.5) with other approaches. 



3.4 The dependence of e'/e on R 6 and R 8 

In addition, it is also of interest to show the values of e'/e for specific values of Rq, Rg and 
A^g-. We do this in table 3. The results in table 3 should be self explanatory. Depending on 
the values of the parameters R e and R s , one finds values of e'/e that are within or outside 
the one standard deviation experimental range given in (1.1). However, the remarkable feature 
of the values in this table is that most of them are within two standard deviations from the 
experimental value. 
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Rq - 


= 1.0 


Rq - 


= 1.2 


Rs = 


= 1.4 




R$ = 0.8 


R 8 = 1.0 


R 8 = 0.8 


R 8 = 1.0 


R 8 = 0.8 


R s = 1.0 


310 


12.2 


10.5 


16.4 


14.8 


20.7 


19.0 


340 


13.3 


11.5 


17.9 


16.1 


22.5 


20.7 


370 


14.6 


12.6 


19.5 


17.6 


24.5 


22.5 



Table 3: The ratio e'/e in units of 10 4 for wit = 165 GeV and various values of ALL Rq, Rg. 

4 Comparison with other Analyses 

Let us now compare our findings of (3.3), (3.4) and (3.5) with other recent estimates of these 
quantities. 

In the approach of refs. [10, 12, 15], the enhancement of the QCD penguin matrix element 
(Qq)o above the factorisation result was attributed to the strong final-state interactions in the 
final state pion system. In [12,15], the final-state interactions were taken into account through 
dispersion relation techniques which resulted in an Omnes type exponential, giving an enhance- 
ment factor 1.55 ± 0.10 which provides full agreement with the experimental result for e'/e. 1 
No such enhancement has been found in the case of (Qs)2, for which a mild suppression factor 
0.92±0.03 was obtained [10]. Furthermore, the results of the Omnes approach were corroborated 
by an explicit computation within chiral perturbation theory at one loop. 

However, the strange mass used in ref. [10] corresponds to (m s + md)(m c ) = 135 MeV, thus 
being 10% higher than the value advocated in the present work. Translating the findings of [10] 
for the case of Omnes resummed final-state interactions into Rq and Rg, we obtain 

Re = 1.05 ± 0.06 , R 8 = 0.84 ± 0.05 , (4.1) 

in reasonable agreement with our values, albeit pointing to a slightly different scenario as to their 
composition into £>-parameter and strange mass. With the present uncertainty in the strange 
quark mass, which is not included in the errors given in eq. (4.1), we are, however, unable to 
decide which is the correct scenario. The central result of (4.1) is also indicated as the full circle 
in figure 1. 

In the very recent article [9], the matrix element of the QCD penguin operator Qq was 
calculated within an approach which was able to take into account the dominant 0(nf/N c ) cor- 
rections to the large- result, where nj is the number of active quark flavours. In addition, the 
authors demonstrated an explicit cancellation of the renormalisation scale dependence between 
1 Further discussion of this approach can also be found in refs. [57,58]. 
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Wilson coefficients and matrix elements. Numerically, in our language, the results of [9] and [53] 
would correspond to 



# 6 = 2.1 ±1.1, i? 8 = 2.2±0.4. (4.2) 

As seen in figure 1, the central result, displayed as the full diamond, is roughly consistent with 
the experimental value for e'/e, but lies far outside the ellipse found here. However, in the 
approach of [9], the dominant uncertainty for Rq is due to the value of the quark condensate, 
and, as already remarked above, both parameters are presently only calculated in the chiral 
limit. Thus, in our opinion, at the moment it is premature to draw any definite conclusions 
from the results of [53] and [9]. 

Similar results of a rather large value of the matrix element (Qe)o were also found earlier in 
the framework of the l/N c expansion at NLO with chiral perturbation theory at the leading order 
[11, 13]. With the value m s (2 GeV) = 105 MeV advocated there, the finding B^ /2) = 2.5 ± 0.4 
leads to Rq = 2.2 ± 0.4, in perfect agreement to the result of eq. (4.2). Together with the result 
Rg = 1.1 ± 0.3 [54], in figure 1 the findings of [11, 13,54] are displayed as the full triangle. 
It should be clear that this set of parameters leads to a value of e'/e much larger than the 
experimental average. Nevertheless, one should note that like in ref. [9], also in [11,13], the 
penguin matrix element (Qe)o was expressed in terms of the quark condensate. It is, however, 
well known that the chiral corrections to the GMOR relation which relates the parametrisations 
(2.10), and the one in terms of the quark condensate, in the strange quark case are large, of 
order 50% [48]. This can easily lead to a factor of two difference for both parametrisations, and 
one has to wait for a complete computation at the next-to-leading order, before it is clear which 
one is closer to the true result. 

To conclude this section, we further elaborate on the strict large-iV c limit and perform the 
following exercise. We set 

B^ /2) = b£ /2) = 1.0 [27,59] and B K = 0.75 [60,61], (4.3) 

In that case one finds [35] 

IniA t = (1.34 ± 0.06) • 10" 4 . (4.4) 

Setting lmX t to this central value, in table 4 we show the results for e'/e as a function of m s {m c ) 
and Ag> . 

We observe that for m s {m c ) m 105 MeV or equivalently m s (2GeV) « 90 MeV and rs 
340 MeV, the results in table 4 are fully consistent with the data on e'/e. The required value of 
m s is smaller, but consistent with (3.2) and in the ballpark of values found in the most recent 
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m s (m c ) = 115 MeV 


m s (m c ) = 105 MeV 


m s (m c ) = 95 MeV 


310 


10.8 


13.4 


16.8 


340 


11.8 


14.6 


18.2 


370 


12.9 


15.9 


19.9 



Table 4: The ratio e'/e in units of 10~ 4 for the strict large-iV c results Bf' ' = Bf l] = 1.0, 
Bk = 0.75 and various values of and m s (m c ). 

lattice simulations with dynamical fermions [44,62,63]. As can be seen from the dashed-dotted 
line in figure 1, the required value for the large- N c condition Rq = R$ is 

R 6 = R 8 = 1.36 ± 0.3 , (4.5) 

which is close to the ellipse found here but outside of it. 



5 Conclusions 

In this letter, we have reanalysed the ratio e'/e in view of the improved experimental data and 
an improved information on the relevant theoretical input parameters. After the uncertainties 
on ImXt have been reduced significantly during the last four years and some consensus has 
been reached on the matrix element of the dominant electroweak penguin operator Q$, the 
main uncertainty in the estimate of e'/e in the SM is the matrix element of the leading QCD 
penguin operator Qq. Assuming that e'/e is fully dominated by the SM contributions, we have 
determined the ratio Rq and the matrix element (Q6)o DR ( m c) from the data on e'/e with the 
results given in (3.4) and (3.5). 

We have also shown that with Bk = 0.75 and m s (m c ) m 105 MeV, the large-A^ values of 
(Qe)o and (Qs}2 with their Wilson coefficients calculated within the SM are consistent with 
the experimental data on e'/e. This should be contrasted with the situation of 1999, when the 
large- N c estimates of e'/e were by a factor of two to three lower than the data. The reasons for 
this change have already been listed in section 1. 

Finally, we have shown in figure 1 the correlation between Rq and Rs as implied by the 
experimental data on e'/e. 

There are three messages from our paper: 

• If indeed Rs = 1.0 ± 0.2 as indicated by several recent estimates, then the data on e'/e 
imply 

Re = 1.23 ±0.16, = 1.23± 0.29 , (Q 6 )o DR (m c ) = - (0.77 ± 0.19) (Q 8 )^ DR (m c ) • 
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(5.1) 



This is in accordance with the results in [10, 12, 15], but as seen in figure 1 it differs from 
the large-iV c approaches in [27-29], [11, 13,54] and [9] with the first one closest to our 
findings. 

• The large-iV c approach of [27-29] can only be made consistent with data provided 

R 6 = R S = 1.36 ± 0.3 (5.2) 

and (Qs}2 DR ( m c) is larger than obtained by most recent approaches discussed above. 
This would require m s {m c ) rj 105 MeV which is still compatible with (3.2). Moreover, 
lower values of m s (m c ) are indicated by the recent lattice simulations with dynamical 
fermions [44,62,63], but the situation certainly requires further study. 

• Higher values of (Q6)o' DR ( ?71 c) and (Qs^^i 771 ^) are obtained in the approach of [9] but 
as the calculations have been done in the chiral limit, it would be premature to draw 
definite conclusions at present. From the point of view of the approach in [27-29] the 
values Rq ~ R$ ~ 2 would correspond to rather small values of m s (m c ) so that a better 

(1/2) (3/2) 

interpretation would be a significant departure of the parameters Bq and B% from 
unity. As the analysis [9] goes beyond the strict large-iV c limit, this is certainly conceivable. 

As shown in figure 1, all these three scenarios are consistent with the data. Which of these pic- 
tures of e' ' je is correct, can only be answered by calculating (Qq}o by means of non-perturbative 
methods that go beyond the large-A^ limit and are reliable. Moreover, the calculation of the 
strange quark mass m s and of {Qs)i should be improved. Such calculations are independent 
of the assumption about the role of new physics in e'/e that we have made in order to extract 
(Qq)o from the data. 

If the values for these matrix elements and i?6,8 will be found one day to lie significantly 
outside the allowed region in figure 1, new physics contributions to e'/e will be required in order 
to fit the experimental data. 
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